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Abstract 

< 

A model for N2 Vegard-Kaplan band (__ 3 __+ — emissions in Martian dayglow has been developed to 

explain the recent observations made by the SPICAM ultraviolet spectrograph aboard Mars Express. Steady 

PL 

state photoelectron fluxes and volume excitation rates have been calculated using the Analytical Yield Spectra 
(AYS) technique. Since inter-state cascading is important for triplet states of N 2 , the population of any given 
level of N2 triplet states is calculated under statistical equilibrium considering direct excitation, cascading, and 
quenching effects. Relative population of all vibrational levels of each triplet state is calculated in the model. 
Line of sight intensities and height-integrated overhead intensities have been calculated for Vegard-Kaplan 
(VK), First Positive (B 3 U g - A 3 E+), Second Positive (C 3 II U - B 3 U g ), and Wu-Benesch (W 3 A U - B 3 U g ) 
bands of N 2 . A reduction in the N 2 density by a factor of 3 in the Mars Thermospheric General Circulation 
-__, Model is required to obtain agreement between calculated limb profiles of VK (0-6) and SPICAM observation. 

Calculations are carried out to asses the impact of model parameters, viz., electron impact cross sections, solar 
EUV flux, and model atmosphere, on the emission intensities. Constraining the N2/CO2 ratio by SPICAM 
observations, we suggest the N2/CO2 ratios to be in the range 1.1 to 1.4% at 120 km, 1.8 to 3.2% at 140 km, 
and 4 to 7% at 170 km. During high solar activity the overhead intensity of N2 VK band emissions would be 
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~2.5 times higher than that during low solar activity. 

1 Introduction and Dalgarno, 1979), who suggested that a high reso- 

. ^ lution UV spectrometer could detect the N2 UV emis- 

^ Emissions from excited states of N 2 have been stud- sions on Mars _ Fox and Da i garno (1979) have pre- 

^ ied extensively in the terrestrial airglow and aurora dicted the intensity of various N 2 triplet state emis- 

(e.g., Sharp, 1971; Conway and Christensen, 1985; sions (Vegard-Kaplan, First positive, Second Positive, 

Meier, 1991; Morrill and Benesch, 1996; Broadfoot W _ B ), along with LBH band of N 2 and First Neg- 

et al, 1997). But the absence of any emission fea- ative band of N +_ In the terrestria i atmosphere emis- 

ture of N 2 during Mariner observations of Mars (Barth sion from Vegard-Kaplan (VK) bands are weak due 

et al, 1971) surprised the planetary scientists who at- to efficient quen ching by atomic oxygen, but C0 2 is 

tributed it to the low fractional abundance by volume not good at quench i n g VK bands (Fox and Dalgarno, 

of molecular nitrogen on Mars {Dalgarno and McEl- mg . Dreyer et ^ m ^ gQ the intensity of these 

roy, 1970). Earlier, N 2 emissions on Mars were pre- band ghould be ap p recia ble in Martian airglow. 
dicted by Fox and co-workers (Fox et al., 1977; Fox 
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Recent observations by SPICAM (Spectroscopy 
for Investigation of Characteristics of the Atmosphere 
of Mars) onboard Mars Express (MEX) have, for the 
first time, observed N2 emissions in the dayglow of 
Mars (Leblanc et al., 2006, 2007). The main emis- 
sions observed are (0, 5), (0, 6), and (0, 7) bands of 
VK system, which originate from triplet A 3 T,+ state 
of excited N2 molecule. The overhead intensity of the 
N2 VK (0, 6) band derived from the intensity observed 
by the SPICAM is found to be ~3 times smaller than 
the intensity calculated by Fox and Dalgarno (1979). 

There have been several measurements of elec- 
tron impact cross sections of triplet states of N2 since 
Fox and Dalgarno (1979) carried out their calcula- 
tions. With new cross sections and updated molec- 
ular parameters (transition probability and Franck- 
Condon factor) a model of N2 dayglow emission on 
Mars is necessary for a better understanding of the 
recent SPICAM observations. In the present work, a 
model has been developed to calculate the N2 day- 
glow emissions on Mars using the Analytical Yield 
Spectra approach. While calculating the emission 
of VK bands of N2, cascading from the higher ly- 
ing states and quenching by atmospheric constituents 
are considered and the population of any given vi- 
brational level of a state is calculated under statisti- 
cal equilibrium. Height-integrated overhead intensi- 
ties are reported for major vibrational bands of N2 
VK, First Positive (B 3 U g - A 3 T,+), Second Positive 
(C 3 U U - B 3 U g ), and Wu-Benesch (W 3 A U - B 3 U g ) 
bands. Limb profiles of VK (0, 5), (0, 6), and (0, 7) 
bands are calculated. The limb profile of VK (0, 6) 
band is compared with that reported by the SPICAM 
observations. The present model has been used re- 
cently to estimate the N2 triplet band intensities on 
the Venus (Bhardwaj and Jain, 2011). 



(X X S+) to the A 3 S+ state is dipole forbidden, so pho- 
toelectron impact is the primary excitation source for 
this state. In addition to the direct excitation from 
the ground state, cascade from higher triplet states 
C, B, W, and B' are also important. All excitations 
of higher triplet states will eventually cascade into 
the ^4 3 S+ state (Cartwright et al., 1971; Cartwright, 
1978). The width and shape of VK bands are quite 
sensitive to the rotational temperature, making them 
useful as a monitor of the neutral temperature of the 
upper atmosphere (Broadfoot et al., 1997). 

All transitions between the triplet states of N2 and 
the ground state are spin forbidden, therefore excita- 
tion of these states is primarily due to the electron 
impact. The higher lying states C, W , and B' pop- 
ulate the B state, which in turn radiates to the A 
state. Inter-system cascading -B 3 II g ^ ^4 3 ^+ and 
B 3 H g f± VF 3 A U is important in populating the B 
state (Cartwright et al, 1971; Cartwright, 1978). 

Direct excitation of the 1/ = vibrational level of 
the A 3 T,+ state by electron impact is extremely small, 
because Frank-Condon factor to the v" = level of 
the ground electronic state, qoo, is only 9.77 x 10 
(Gilmore et al., 1992; Piper, 1993). Contributions 
to u' = level of A state come from the higher 
states cascading. We have also included E — » B, 
E —7- C ', E —> A, B —7- W, and reverse first pos- 
itive A — > B cascading in our calculation. The ef- 
fect of reverse first positive transition is important 
in populating the lower vibrational levels of B state, 
which in turn populate the lower vibrational levels 
of the A state (Sharp, 1971; Cartwright et al, 1971; 
Cartwright, 1978). Thus, to calculate the production 
rate of any vibrational level of triplet state of N2 , one 
must take into account direct excitation as well as 
inter-state cascading effects. 



2 Vegard-Kaplan Band (^ 3 S+ ->• 

Figure 1 shows schematic diagram of N2 triplet states 
energy level with excitation and subsequent cascad- 
ing processes. The transition from the ground state 



3 Model Input Parameters 

The model atmosphere considering five gases (CO2, 
CO, N2, O, and O2) is taken from the Mars Ther- 
mospheric General Circulation Model (MTGCM) of 
Bougher et al. (1990, 1999, 2000) for a solar longi- 
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tude of 180°, latitude of 47.5°N, and at 1200 LT; and 
is same as used in the study of Shematovich et al. 
(2008). The EUVAC model of Richards et al. (1994) 
has been used to calculate the 37-bin solar EUV flux 
for the day of observation, which is based on the F10.7 
and F10.7A (81-day average) solar index. The F10.7 
flux as seen by Mars (by accounting for the Mars-Sun- 
Earth angle) is used to derive the 37-bin solar EUV 
flux. The EUVAC solar spectrum thus obtained is 
then scaled for the heliocentric distance of Mars for 
the day, considered in the present study. To assess 
the impact of solar EUV flux on model calculations, 
we have also used SOLAR2000 v.2.36 (S2K) model of 
Tobiska (2004). 

Photoionization and photoabsorption cross sec- 
tions for the gases considered in the present study are 
taken from Schunk and Nagy (2000). The branching 
ratios for excited states of CO^~, CO + , Ng , + , and 

have been taken from Avakyan et al. (1998). For 
calculating the intensity of a specific band v' — v" , 
Franck-Condon factors and transition probabilities 
are required. For N2 these are taken from Gilmore 
et al. (1992). Electron impact cross sections for N2 
triplet excited states (A, B, C, W, B' , and E) were 
measured by Cartwright et al. (1977) up to 50 eV. 
These cross sections were renormalized later by Tra- 
jmar et al. (1983) with the use of improved data on 
elastic cross sections. More recently, N2 triplet state 
cross sections have been measured by Campbell et al. 
(2001) and Johnson et al. (2005). Itikawa (2006) re- 
viewed the cross sections of the N2 triplet excited 
states and recommended the best values determined 
by Brunger et al. (2003). We have taken the N2 triplet 
states cross sections from Itikawa (2006), which have 
been fitted analytically using equation (cf. Jackman 
et al., 1977; Bhardwaj and Jain, 2009) 
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where qo = A-koqR 2 and has the value 6.512 x 10~ 14 
eV 2 cm 2 . Table 1 shows the corresponding parame- 
ters. Fig. 2 shows the fitted cross sections of the N2 



triplet A, B, C, and W states along with the recom- 
mended cross sections of Itikawa (2006). For other 
gases electron impact cross sections have been taken 
from Jackman et al. (1977), except for CO2, which 
are from Bhardwaj and Jain (2009). 

We have run our model for the Mars Express ob- 
servation on 16 Dec. 2004 (Sun-Mars distance = 1.59 
AU, and F10.7 at Mars = 35.6), taking solar zenith 
angle as 45°, solar EUV flux from the EUVAC model, 
and MTGCM model atmosphere. Hereafter we refer 
it as the "standard case" . We have also studied the ef- 
fects of various input parameters (like solar EUV flux, 
N2 triplet state cross sections, model atmosphere, so- 
lar cycle) on the emission intensity, which are dis- 
cussed in Section 6. 



4 Model Calculation 

4.1 Photoelectron Production Rate 

Primary photoelectron production rate is calculated 
using 



1 i,A 



(he 



I(Z, A) = I(oo, A) exp 



sec(x)XX( A ) 



- E — Wji 

(2) 
ni{Z )dZ 

(3) 

where af and af(j,\) are the total photoabsorption 
cross section and the photoionization cross section of 
the jth ion state of the constituent I at wavelength A, 
respectively; I(oo, A) is the unattenuated solar flux at 
wavelength A, n\ is the neutral density of constituent 
I at altitude Z; sec(x) is the Chapman function, % 
is the solar zenith angle (SZA); 5(hc/X — E — Wji) 
is the delta function, in which hc/X is the incident 
photon energy, Wji is the ionization potential of jth 
ion state of the Zth constituent, and E is the energy 
of ejected electron. We have used sec(x) in place of 
ch(x), which is valid for % values upto 80°. Figure 3 
shows the primary photoelectron energy spectrum at 
three different altitudes. There is a sharp peak at 27 
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eV due to the ionization of CO2 in the ground state 
by the He II solar Lyman a line at 303.78 A. The 
peaks at 21 and 23 eV are due to ionization of CO2 
in the A 2 II M and B 2 £+ states of CO^", respectively, 
by the 303.78 A solar photons. The individual peaks 
structure shown in the figure are different from that of 
Mantas and Hanson (1979), which is due to revisions 
in the branching ratios (Avakyan et al, 1998) used in 
the present study. 



U(E, Eo) = A lE ° k + ME^/e^n + 

(5) 

Here E k = £ /1000, e = E/I (I is the lowest ion- 
ization threshold), and x = (E — B<i)jB\. A\ = 
0.027, A 2 = 1.20, t = 0, r = 0, s = -0.0536, 
B = 10.095, Bi = 5.5, and B 2 = 0.9 are the best 
fit parameters. 

For other gases, viz., O2, N2, O, and CO, we have 
used the AYS given in Singhal et al. (1980) 



4.2 Photoelectron Flux 

To calculate the photoelectron flux we have adopted 
the Analytical Yield Spectra (AYS) technique (cf. 
Singhal and Haider, 1984; Bhardwaj and Singhal, 
1990; Bhardwaj et al., 1990, 1996; Singhal and Bhard- 
waj, 1991; Bhardwaj, 1999, 2003; Bhardwaj and 
Michael, 1999a,b). The AYS is the analytical rep- 
resentation of numerical yield spectra obtained us- 
ing the Monte Carlo model (cf. Singhal et al., 1980; 
Bhardwaj and Michael, 1999a, b; Bhardwaj and Jain, 
2009). Recently, the AYS model for electron degra- 
dation in CO2 has been developed by Bhardwaj and 
Jain (2009). Further details of the AYS technique are 
given in Bhardwaj and Michael (1999a), Bhardwaj and 
Jain (2009), and references therein. Using AYS the 
photoelectron flux has been calculated as (e.g. Singhal 
and Haider, 1984; Bhardwaj and Michael, 1999b) 

«Z,W m ^t**U (4) 
Jw kl 2_^ni(Z)a lT (E) 

1 

where <jit{E) is the total inelastic cross section for 
the lib. gas, n\ is its density, and U(E, Eq) is the two- 
dimensional AYS, which embodies the non-spatial in- 
formation of degradation process. It represents the 
equilibrium number of electrons per unit energy at an 
energy E resulting from the local energy degradation 
of an incident electron of energy Eq. For the CO2 gas 
it is given as (Bhardwaj and Jain, 2009) 



U(E,E ) = C + C^Ek + K)/[(E - M) 2 + L 2 }. (6) 

Here Co, C\, K, M, and L are the fitted parameters 
which are independent of the energy, and whose values 
are given by Singhal et al. (1980). 

The calculated photoelectron flux at 130 km al- 
titude is shown in Figure 4 for the standard case as 
well as for conditions similar to those of Viking 1 (see 
Section 5.2). The photoelectron flux calculated by Si- 
mon et al. (2009) and Fox and Dalgarno (1979) are 
also shown in Figure 4 at same altitude. Overall im- 
portant peak structures are similar in all the three 
calculated fluxes, e.g., the peak at 27 eV and broad 
peak at 21-23 eV. A sharp dip at around 3 eV is promi- 
nent in all three photoelectron fluxes, which is due to 
large vibrational cross sections at 3.8 eV for electron 
impact on CO2. The calculated fluxes decrease expo- 
nentially with increasing energy. The sudden decrease 
in the photoelectron flux at higher energies is due to 
the presence of these features in the primary photo- 
electron energy spectrum (cf. Figure 3). 

5 Results and discussion 

5.1 Volume excitation rates 

We have calculated volume excitation rate Vu(Z,E) 
for the ith state of the Ith gas at altitude Z and en- 
ergy E using the equation (Singhal and Bhardwaj, 
1991; Bhardwaj, 1999, 2003; Bhardwaj and Michael, 
1999b) 
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V u (Z,E)=m(Z) (j>{Z,E)a il {E)dE, (7) 

where ni(Z) is the density of the Ith. gas at altitude 
Z and uu{E) is the electron impact cross section for 
the ith state of the lib. gas, for which the threshold is 
E t h- Figure 5 (upper panel) shows the volume exci- 
tation rates of the N2 triplet states (^4, B, C, W, B' , 
and E) excited by photoelectron impact. The altitude 
of peak production for all states is ~126 km for the 
standard case. The volume excitation rate of N2(A) 
state calculated using the S2K solar flux model is also 
shown in the upper panel of Figure 5. The peak of ex- 
citation rate occurs at the same altitude for both solar 
EUV flux models but the magnitude of excitation rate 
is slightly higher when the S2K model is used. More 
discussion about the effect of solar EUV flux model 
on emission intensities is given in Section 6.2. 

To calculate the contribution of cascading from 
higher triplet states and interstate cascading between 
different states, we solve the equations for statistical 
equilibrium based on the formulation of Cartwright 
(1978) and assumed that only excitation from the low- 
est vibrational level of the electronic ground state is 
important. At a specified altitude, for a vibrational 
level v of a state a, the population is determined using 
statistical equilibrium 

FV + ^^Af>? = {if; + ^^A^K (8) 

/3 s 7 r 

where 

V a electron impact volume excitation rate 

(cnr 3 s _1 ) of state a; 
qo u Franck-Condon factor for the excitation 

from ground level to v level of state a; 
Asu transition probability (s™ 1 ) from state 

f3(s) to a(y); 
Kq V total electronic quenching frequency 

(s _1 ) of level v of state a by the 

all gases defined as: Yl^qn) v x n l'i where, 



^q(i)u * s ^ ne quenching rate coefficient 

of level v of a by gas I of density n\ ; 
Ayr transition from level v of state a to 

vibrational level r of state 7; 
n density (cm -3 ); 

a, (3, 7 electronic states; 
s, r source and sink vibrational levels, 

respectively. 

While calculating the cascading from C state, we 
have taken predissociation also into account. The C 
state predissociates approximately half the time (this 
is an average value for all vibrational levels of the C 
state; and 1 levels do not predissociate at all) (cf. 
Daniell and Strickland, 1986). In the terrestrial ther- 
mosphere, the N2(A) state is effectively quenched by 
atomic oxygen. In the case of Mars the main con- 
stituent CO2 does not quench N2(A) level that effi- 
ciently, but still there will be some collisional deactiva- 
tion by other atmospheric constituents of Mars. The 
electronic quenching rates for vibrational levels of N2 
triplet states by O, O2, and N2 are adopted from Mor- 
rill and Benesch (1996) and Cartwright (1978) and by 
CO2 and CO are taken from Dreyer et al. (1974). 

Figure 6 shows the population of different vibra- 
tional levels of triplet states of N2 relative to the 
ground state at 130 km. The relative population of 
N2(A) at 110 km is also shown in the figure. Our 
calculated relative vibrational populations agree well 
with the earlier calculations (Morrill and Benesch, 
1996; Cartwright, 1978). To show the effect of quench- 
ing the relative vibrational populations of N2 (A) state 
calculated without quenching at 110 and 130 km are 
also shown in Figure 6. The quenching does affect the 
vibrational population of N2(A) state mainly for vi- 
brational levels between 5 and 10 at lower altitudes 
(<130 km), as the altitude increases the effect of 
quenching decreases. Figure 7 shows the steady state 
fractional population altitude profiles of a few vibra- 
tional levels of A state and u' = level of B, C, W, 
and B' excited states of N2. 

After calculating the steady state density of differ- 
ent vibrational levels of excited states of N2, the vol- 
ume emission rate V^„ of a vibration band v' — > v" 
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can be obtained using 

Vf u „ = n a u , x A a J u „ (cm^s- 1 ) (9) 

where n®, is the density of vibrational level v' of state 
a, and A^ v „ is the transition probability (s -1 ) for the 
transition from the v' level of the a state to the v" 
level of the (3 state. Figure 5 (bottom panel) shows 
the volume emission rates for the VK (0, 4), (0, 5), 
(0, 6), and (0, 7) bands. We have integrated the vol- 
ume emission rates over the altitudes 80 — 400 km to 
obtain the overhead intensity for various VK bands of 
N2, which are tabulated in Table 2 (standard case). 
Table 3 shows the calculated height-integrated over- 
head intensities for a few of the prominent bands of 
First Positive (B -> A), Second Positive (C B), 
and Wu-Benesch (W — > B) emissions. 

5.2 Line of sight intensity 

For comparison of the calculated intensity with SPI- 
CAM observation we have integrated the calculated 
emission rate along the line of sight and expressed the 
results in kR (1 Rayleigh = 10 6 photon cm~ 2 s _1 ) 

I = J v(r)dr, (10) 

where V(r) is the volume emission rate (in cm~ 3 s _1 ) 
for a particular emission, calculated using equation (9) 
and r is abscissa along the horizontal line of sight. 
The upper limit of the atmosphere in our model is 
taken as 400 km. While calculating limb intensity we 
assume that the emission rate is constant along lo- 
cal longitude/latitude. For the emissions considered 
in the present study, the effect of absorption in the 
atmosphere is found to be negligible. As mentioned 
earlier (cf. Leblanc et al., 2007), the main N2 emis- 
sion features observed by SPICAM are (0, 5) and (0, 
6) transitions of the Vegard-Kaplan (A 3 £+ — X 1 ^^) 
band. Leblanc et al. (2006) also reported the detec- 
tion of VK (0, 7) band, but it was characterized by a 
large uncertainty because it falls between two intense 
emissions at 289 nm and 297.2 nm of CO^ UV doublet 
and oxygen line emission, respectively. Otherwise, as 



shown in Table 2, VK(0, 7) band would have been 
more intense than the (0, 5) band. The ratio between 
calculated intensity of the VK (0, 6) and (0, 5) bands 
is 1.3, which is in good agreement with the results of 
Leblanc et al. (2007) and Fox and Dalgarno (1979). 

Figure 8 shows the limb profiles of the VK (0, 6) 
band at different solar zenith angles along with the 
SPICAM observed profiles averaged over the solar lon- 
gitude L s 100°-171° and SZA 8°-36° and 36°-64°, 
taken from Leblanc et al. (2007). The effect of SZA 
on the calculated profiles is clearly visible in Figure 8; 
the peak of the altitude profile rises while the inten- 
sity decreases with increasing SZA. The limb profiles 
of the VK (0, 5) and (0, 6) bands at SZA=45° are also 
plotted in Figure 8. For the standard case (SZA=45°), 
the peak intensities of the VK (0, 5), (0, 6), and (0, 
7) bands are ~0.9, 1.1, and 1 kR, respectively, at 120 
km. For SZA values of 20° and 60°, the N 2 VK (0, 
6) band peaks at 118 and 124 km with a value of 1.4 
and 0.9 kR, respectively. 

The shape of calculated and observed limb inten- 
sities are in agreement with each other but the magni- 
tude of calculated intensities are larger by a factor of 
~3 at SZA = 20°. This difference could be due to the 
larger abundance of N2 in the model atmosphere used 
in the present study. Other factors can also affect the 
calculated intensities, but their combined uncertain- 
ties also cannot account for the difference by a factor 
of 3 in the calculated and observed intensities (effect 
of other input parameters, viz., electron impact cross 
section, and solar EUV flux model is described in Sec- 
tions 6.1 and 6.2, respectively). Figure 8 also shows 
the computed limb intensity of the VK (0, 6) emis- 
sion at SZA 20°, 45°, and 60° obtained after reducing 
the density of N2 by a factor of 3, which compares 
favourably in both shape and magnitude with the ob- 
served emission. The N2/CO2 ratio, after reducing 
the N2 density by a factor of 3 is to 0.9, 2.1, and 7.1% 
at altitudes of 120, 140, and 170 km, respectively. The 
calculated overhead intensities of VK bands after re- 
ducing N2 density by a factor of 3 (for the standard 
case) are depicted in column 3 of Table 2. It may how- 
ever be noted that the observed limb profiles (Leblanc 
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et al., 2007) are averaged over several days of observa- 
tion (L s =101°-171°) and range of SZA values, while 
the model profile is for a single day (16 Dec. 2004) at 
L s = 130° and SZA = 20°. 

We have also calculated the nadir intensity for the 
condition similar to that of Viking landing (Sun-Mars 
distance = 1.65 AU and F10.7 = 68). The model at- 
mosphere was taken from Fox (2004) for the low solar 
activity condition and a SZA of 45°. For the VK (0, 
6) band our calculated intensity is 26 R, which is con- 
sistent with results (20 R) of Fox and Dalgarno (1979) 
for the similar condition. The minor difference may be 
due to the updated cross sections and transition prob- 
abilities. For the condition similar to that of Viking, 
Leblanc et al. (2007) have measured an intensity of 
~180 R for the VK (0, 6) band, which corresponds 
to a nadir intensity of ~6 R. The measured value is 
about 4 times smaller than our calculated intensity. 
Such a difference by a factor of 4 between observed 
and calculated intensities might be due to the higher 
density of N2 taken in our model atmosphere. Leblanc 
et al. (2007) mentioned that difference by factor of 3 
between the estimated and nadir intensity calculated 
by Fox and Dalgarno (1979) could have been due to 
the larger N2/CO2 ratio in the model atmosphere of 
Fox and Dalgarno (1979). Leblanc et al. (2007) have 
suggested that ratio of the integrated column densi- 
ties of N2 and CO2 between 120 and 170 km, that is 
the mixing ratio between N2 and CO2 for a uniformly 
mixed atmosphere, would be 0.9% for an overhead in- 
tensity of 6 R. For the same altitude range, the ratio 
of N2/CO2 density is 3.5% and 3.7% in model atmo- 
sphere used in the work of Fox (2004) and Bougher's 
MTGCM, respectively, which is a factor of 4 higher 
than that suggested by Leblanc et al. (2007). 

To summarize, the above results indicate that the 
N2 density in the MTGCM atmosphere, as well as in 
the model atmosphere of Fox (2004), has to be re- 
duced by a factor of ~3 to obtain agreement between 
the SPICAM observation and the calculated intensity. 



5.3 Variation with Solar Zenith Angle and 
Solar 10.7 flux 

Figure 9 shows the variation of the VK (0, 6) band in- 
tensity, averaged between 120 and 170 km, with SZA 
and its comparison with SPICAM observations. Cal- 
culated intensities are for standard case obtained after 
reducing the N2 density profile in the MTGCM atmo- 
sphere by a factor of 3 (see discussion in the previous 
section). Model intensity shows a cosine SZA depen- 
dence, with larger attenuation of solar EUV flux at 
higher SZA, resulting in decrease in the intensity at 
higher SZA. Calculated intensities are in agreement 
with the observed values, within observational and 
model uncertainties. 

Another important model parameter, which af- 
fects the emission intensities is the solar EUV flux, 
whose variation is assumed to be given by the F10.7 
index. Solar EUV flux has been calculated using the 
F10.7 flux for the day of observation and scaled to 
the Mars according to its heliocentric distance. For 
the observations reported by Leblanc et al. (2007) the 
solar longitude of Mars varied between 101° and 171°, 
which corresponds to change in the heliocentric dis- 
tance of Mars form 1.64 to 1.49 AU. Figure 10 shows 
the variation of VK (0, 6) band intensity with respect 
to the F10.7 solar index at Mars. Calculations are 
made for the standard case with the N2 density in the 
MTGCM model reduced by a factor 3. Model calcu- 
lated intensities are consistent with the observed val- 
ues within the uncertainties of observation and model. 



6 Effect of various model parame- 
ters on Intensity 

To evaluate the effect of various model input param- 
eters, such as solar flux, cross sections, and model 
atmosphere, on the VK band emissions, we have con- 
ducted a series of test studies by changing one param- 
eter at a time and compare the results with those of 
the standard case. The results are presented in Ta- 
ble 2 and discussed below. 
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6.1 Electron impact cross sections for the 
triplet states 

Since electron impact on N2 is the source of excita- 
tion of forbidden triplet states of N2, any change in 
electron impact cross sections will directly affect the 
VK band emission intensities. Various measurements 
of the N2 triplet state cross sections were discussed 
in Section 3. In the standard case we have taken the 
recommended cross sections of Itikawa (2006), which 
are fitted using the semiempirical relation given in 
equation (1) (cf. Table 1 and Figure 2). Instead of 
analytically fitted cross sections, if the triplet state 
cross sections of Itikawa (2006) are used in the model, 
the calculated triplet band intensities differ from the 
standard case by less than 10%. 

Itikawa's recommended cross sections are based 
on the best values determined by the Brunger et al. 
(2003). For the triplet states cross section, Brunger 
et al. (2003) have estimated the uncertainty of the 
recommended cross sections as ±35% (±40% at en- 
ergies below 15 eV) for A 3 £+, ±35% for B 3 U g and 
W 3 A U , ±40% for B' 3 Z~, ±30% for C 3 U U and ±40% 
for E 3 T,^ state. The integral cross sections (ICS) of 
Johnson et al. (2005) are derived from the differential 
cross sections (DCS) of Khakoo et al. (2005). Johnson 
et al. (2005) have given the ICS at 8 energies between 
10 and 100 eV, with uncertainty for all states cross 
sections varying between ±20% to ±22%; at a few 
energy points it is as high as ±35%. 

To evaluate the effect of electron impact cross sec- 
tions on the VK band emissions we have taken two sets 
of cross sections; one from Cartwright et al. (1977), 
which were renormalized by Trajmar et al. (1983), 
and second from the recent cross sections given by 
Johnson et al. (2005). The resulting VK band inten- 
sities are shown in Table 2. The intensities calculated 
with the cross sections of Trajmar et al. (1983) are 
almost the same as in the standard case. However, 
when the cross sections of Johnson et al. (2005) are 
used, the VK band intensities are reduces by 45%, 
compared to the intensities computed for the standard 
case, which is due to smaller cross sections of John- 



son et al. (2005). The effect of the smaller triplet 
state cross sections of Johnson et al. (2005) is also 
seen on the limb intensities shown in the Figure 11 
where a reduction in N2 density by a factor 2 is suf- 
ficient to fit the SPICAM observed profile. Thus, the 
electron impact triplet state excitation cross sections 
of N2 also help in constraining the N2 density in the 
model atmosphere. 

6.2 Input solar EUV flux model 

SOLAR2000 model of Tobiska (2004) and EUVAC 
model of Richards et al. (1994) are the two widely 
used solar flux models in the aeronomical calculations. 
In the standard case we have used EUVAC model. To 
see the effect of input solar flux on the VK emissions, 
we conducted a test study by taking the solar EUV 
flux from SOLAR2000 v.2.36 (S2K) model of Tobiska 
(2004) at 37 wavelength bins; the other input param- 
eters remain the same as in the standard case. The 
calculated integrated overhead intensities are shown 
in the Table 2. The calculated intensities of VK bands 
using S2K model are ~15% larger than those calcu- 
lated by using the EUVAC model. This results in the 
requirement of a larger reduction in the N2 density, 
that is, a factor of 3.4 compared to 3.0 for the stan- 
dard case to fit the observed limb profile of the VK 
(0, 6) band. 

6.3 Model atmosphere 

The importance of model atmosphere on the calcu- 
lated intensities has been demonstrated in Section 5.2. 
The N2/CO2 ratio, which describes the abundance 
of molecular nitrogen in the atmosphere of Mars, is 
different in different model atmospheres. For the 
present study we have taken the atmosphere from 
Bougher's MTGCM {Bougher et al, 1990, 1999, 2000) 
as used in study of Shematovich et al. (2008) where 
the N2/CO2 ratio is 2.8, 6.4 and 21% at 120, 140 and 
170 km, respectively. Leblanc et al. (2007) suggested 
that N2/CO2 ratio is higher in the model atmosphere 
used by Fox and Dalgarno (1979). The recent models 
of Krasnopolsky (2002) are characterized by smaller 
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abundances of N2 than that of Fox and Dalgarno 
(1979). The N 2 /C0 2 ratios are 2.6, 3.8, and 8.6% at 
120, 140 and 170 km, respectively in Krasnopolsky' s 
model. 

We have used the model atmospheres of 
Krasnopolsky (2002) and Fox (2004) to study the ef- 
fect of model atmosphere on the VK emission inten- 
sities. Figure 11 shows the calculated limb intensity 
of the VK (0, 6) band for both model atmospheres at 
SZA 20° (all other conditions are similar to the stan- 
dard case). The emission peaks at ~116 km in the 
case of Krasnopolsky (2002), which is almost similar 
to standard case (~118 km). But the emission peaks 
at higher altitude (~123 km) when the model atmo- 
sphere of Fox (2004) is used, which is due to higher 
CO2 abundance in her model. The intensities calcu- 
lated using both models are found to be larger than 
the observed values. To fit the observed limb profile, 
the N2 density in the Krasnopolsky (2002) model has 
to be reduced by a factor of 2.1, the N2/CO2 ratios 
thus become 1.3, 1.8, and 4.4% at 120, 140, and 170 
km, respectively. In the case of Fox (2004) model 
atmosphere, the required decrease in N2 density is a 
factor of 2.5, which corresponds to the N2/CO2 ra- 
tios of 1.1, 1.9, and 5.3% at 120, 140, and 170 km, 
respectively. 

6.4 Solar Cycle 

The solar cycle is approaching higher solar activity, 
and MEX is currently orbiting Mars. We therefore 
hope to observe the effects of higher solar activity on 
the Martian dayglow emissions. Using the EUVAC 
model we have calculated the various N2 triplet band 
emissions for high solar activity conditions similar to 
that of Mariner 6 and 7 flybys when the F10.7 index 
was ~ 190 at 1 AU. The model atmosphere for so- 
lar maximum conditions was taken from Fox (2004); 
other model parameters are same as in the standard 
case. The calculated height-integrated overhead in- 
tensities for the VK bands are presented in the Ta- 
ble 2, and those for the other triplet bands in Table 3. 
The calculated solar maximum intensities are larger 
by a factor of ~1.5 than those of the standard case, 



and ~2.5 times larger than those for Viking condi- 
tions. As mentioned in section 5.2, the calculated and 
observed limb profiles are consistent with each other 
when the N2 density in the atmosphere is reduced by 
a factor of 3. If a similar situation prevails during high 
solar activity conditions, then the calculated intensity 
of N2 VK band system would be smaller by a factor 
of 2 to 3. 

7 Summary 

We have presented models for the intensities of the N2 
triplet band systems in the Martian dayglow. We have 
used the analytical yield spectra technique to calcu- 
late the steady state photoelectron flux, which in turn 
is used to calculate volume excitation rates of N2 VK 
bands and other triplet states. The populations of 
various vibrational levels of the triplet states of N2 
have been calculated considering direct excitation as 
well as cascading from higher triplet states in statisti- 
cal equilibrium conditions. Using calculated emission 
rates the limb profiles of the VK (0, 5), (0, 6), and (0, 
7) bands have been calculated and compared with the 
SPICAM observed limb profile reported by Leblanc 
et al. (2007). The observed and calculated limb pro- 
files of the VK (0, 6) band are in good agreement 
when the N2 density is reduced by a factor of 3 from 
those given by the MTGCM model of Bougher et al. 
(1990, 1999, 2000). Overhead intensities of prominent 
transitions in VK, First Positive, Second Positive, and 
W —7- B bands have been calculated. 

The effect of important model parameters, viz., 
electron impact N2 triplet state excitation cross sec- 
tions, solar flux, solar activity, and model atmo- 
sphere, on emissions have been studied. Changes in 
cross sections of N2 triplet states can alter the cal- 
culated intensity by a factor of ~2. On the other 
hand, the calculated intensities are ~15% larger when 
the SOLAR2000 v.2.36 solar EUV flux model of To- 
biska (2004) is used instead of the EUVAC model of 
Richards et al. (1994). During high solar activity, 
when the F10.7 is similar to those at the times of 
the Mariner 6 and 7 flybys, the calculated intensities 
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are about a factor of 2.5 larger than those calculated 
for the low solar activity conditions of the Viking mis- 
sion. On using the model atmospheres of Fox (2004) 
and Krasnopolsky (2002), a decrease in N2 density in 
their atmospheric model by a factor of 2.5 and 2.1, re- 
spectively, is required to reconcile the calculated VK 
(0, 6) band limb profile with the observed profile. 

The most important parameter that governs the 
limb intensity of VK band is the N2/CO2 ratio. Con- 
straining the N2/CO2 ratio by SPICAM observations, 
for different cases of model input parameters, we sug- 
gest that the N2/CO2 ratio would be in the range of 
1.1 to 1.4% at 120 km, 1.8 to 3.2% at 140 km, and 4 to 
7% at 170 km. Our study suggests that most of the at- 
mospheric models have N2 abundances that are larger 
than our derived values by factors of 2 to 4. Clearly 
there is a need for improved understanding of the Mar- 
tian atmosphere, and the SPICAM observations help 
to constrain the N2 relative abundances. A decrease 
in the N2 densities in the atmospheric models, as sug- 
gested by our calculations, would affect the chemistry 
and other aeronomical processes in the Martian upper 
atmosphere and ionosphere. 
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Table 1: Fitting parameters (equation 1) for N2 triplet state cross sections. 



_, N2 states 

Parameter 





A 3 £+ 


B 3 n, 


c 3 n u 


W 3 A U 


B' 3 £- 


E 3 £+ 


Th* 


6.17 


7.35 


11.03 


7.36 


8.16 


11.9 


O' 


1 nn 


^ nn 

o.uu 




1 ^n 


1 7n 

r . / u 


1 70 


p 


1.55 


2.33 


1.00 


2.30 


1.50 


3.00 


n 


2.13 


2.50 


2.70 


2.60 


2.12 


3.00 


F 


0.20 


0.178 


0.248 


0.378 


0.08 


0.03 


W 


6.99 


7.50 


11.05 


8.50 


8.99 


12.0 


* Threshold in eV. 
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Table 2: N2 Vegard-Kaplan Band (A 3 T,^ — > X 1 ^^) height-integrated overhead intensity for different cases. 



Band Overhead Intensity (R) 



v> - v" 


Origin 

(A) 


Std.* 
case 


/3.0 


Viking 
Cond. 


Cross section 
CS-AT CS-B* 


Flux 

S2K§ 


Max. 1 


0-2 


2216 


1.5 


0.5 


0.9 


1 
1 


1 A 
1 .^t 


1.7 


2.3 


0-3 


2334 


7.2 


2.5 


4.4 


K 
O 


O.o 


8.3 


10.9 


0-4 


2463 


19.4 


6.8 


11.7 


lo. O 


lo.o 


22.2 


29.3 


0-5 


2605 


34.3 


12.1 


20.7 


ZO.O 


TO A 


39.4 


51.8 


0-6 


2762 


43.7 


15.4 


26.3 


OU 




50.1 


66.0 


0-7 


2937 


41.5 


14.6 


25.0 


28.5 


39.2 


47.6 


62.7 


0-8 


3133 


30.7 


10.8 


18.5 


21 


29 


35.2 


46.4 


0-9 


3354 


18 


6.3 


10.8 


12.3 


17 


20.6 


27.0 


1-8 


2998 


25.9 


9.1 


15.5 


18.8 


25.3 


29.6 


38.4 


1-9 


3200 


38 


13.4 


22.8 


27.8 


37.4 


43.7 


56.6 


1-10 


3427 


35.9 


12.7 


21.5 


26 


35.1 


41 


53.2 


1-11 


3685 


24 


8.5 


14.4 


17.5 


23.5 


27.5 


35.6 


2-10 


3270 


12 


4.2 


7.1 


8.8 


11.9 


13.6 


17.3 


2-11 


3503 


24.9 


8.8 


14.8 


18.5 


24.9 


28.5 


36.3 


2-12 


3769 


26.9 


9.5 


16.0 


20 


26.9 


30.8 


39.3 


2-13 


4074 


18.9 


6.7 


11.3 


14 


18.9 


21.7 


27.6 


3-13 


3857 


16.3 


5.7 


9.7 


12.3 


16.5 


18.7 


24.0 


3-14 


4171 


18.1 


6.4 


10.8 


13.7 


18.3 


20.7 


26.7 



* Standard case. See text for details 
^Cross sections taken from Johnson et al. (2005). 
■'•Cross sections taken from Trajmar et al. (1983). 
§ SOLAR2000 model of Tobiska (2004). 

1 Solar maximum flux for condition similar to Mariner 6 flyby (F10.7 ~ 190). 
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Table 3: Calculated height-integrated overhead intensity of N2 triplet emissions. 



Band Band Origin Intensity (R) 
iv'-v") A Std.* Max. 1 " 



First Positive B 3 U g 






0-0 


10469 


fin q 




0-1 


12317 


39 7 


51 9 


0-2 


14895 


Q R 


1 5 ^! 
10.0 


1-0 


8883 




1 4Q 7 


1-2 


11878 


1 7 R 
± 1 .0 


97 8 


1-3 


14201 


1 ^ 
±0 . \j 


91 9 


2-0 


7732 


4fi 9 


79 Q 


9 1 


RfiQ^ 


uu.u 


Q4 9 


2-2 


9905 


11.2 


17.5 


3-1 


7606 


64.3 


100 


3-2 


8516 


1 fi 8 
lO.o 


ZD. 1 


3-3 


9648 


20.3 


31.5 


4-1 


6772 


19.6 


30.5 


4-2 


7484 


49.2 


76.4 


4-4 


9404 


14.8 


22.9 


5-2 


6689 


22.1 


34.4 


5-3 


7368 


26.2 


40.7 


6-3 


6608 


18 


28 


7-4 


6530 


12 


18.5 


Second Positive C 3 ~tt 


u- B% 


7 


0-0 


3370 


32.1 


51 


0-1 


3576 


21.7 


34.5 


0-2 


3804 


8.7 


13.9 


1-0 


3158 


8.3 


13.2 


Wu-Benesch (W 3 A U - 






2-0 


33206 


4 


6.4 


3-0 


22505 


3.4 


5.4 


3-1 


36522 


3 


4.7 


4-1 


24124 


5.1 


8 


5-1 


18090 


4.7 


7.4 


5-2 


25962 


4.5 


7.0 


6-2 


19193 


5.8 


9.0 


7-2 


15281 


4.7 


7.3 


7-3 


20421 


5.1 


8.0 


8-3 


16112 


5.1 


8.0 


9-4 


17024 


4.4 


6.9 



*Standared Case. See text for details 



t Solar maximum flux for condition similar to Mariner 6 flyby. 
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Figure 1: Energy level diagram for the excitation of N2 triplet states and subsequent inter-state cascading 
processes. Solid arrows show the excitation from ground state to higher states, and dashed arrows represent 
the transitions between different states (HK: Herman-Kaplan; IP: First Positive; RIP: Reverse First Positive; 
2P: Second Positive; VK: Vegard-Kaplan band system). Excitation thresholds for all the triplet states are given 
Table 1. 
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Figure 2: The triple states cross sections due to electron impact on N2. Symbols represent the values of Itikawa 
(2006) and the solid curve represents the analytical fits using equation 1. Cross sections of B and W have been 
plotted after multiplying by a factor of 10 and 5, respectively. 
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Figure 3: Primary photoelectron energy distribution at three different altitudes for the standard case. 
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Figure 4: Model steady-state photoelectron flux calculated at 130 km for standard case and for Viking condition. 
Flux calculated by Simon et al. (2009) and Fox and Dalgarno (1979) at 130 km are also shown for comparison. 
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Volume Emission Rate (cm" s" ) 

Figure 5: (Upper panel) The volume excitation rates of various triplet states of N2 by direct electron impact 
excitation for the standard case. Dashed curve shows the excitation rate of A state calculated using S2K model. 
The excitation rate of the E state has been multiplied by a factor of 10. (Bottom panel) The volume emission 
rates of the VK (0, 4), (0, 5), (0, 6), and (0, 7) bands. 
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Figure 6: The relative populations of vibrational levels of different triplet states of N2 with respect to the N2 
density at 130 km. Dashed line with triangle shows the relative vibrational populations of A at 110 and 130 
km, respectively, without considering the quenching. 




Relative population 



Figure 7: Altitude profiles of the relative populations of selected vibrational levels of the N2(A) state, and 
level of the B, B', C, and W states with respect to those of N2(X). Population of B,B', and C have been 
plotted after multiplying by a factor of 10 4 , 10 5 , and 10 6 , respectively. 
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Figure 8: Calculated limb intensity of the N2 VK (0, 6) band at different solar zenith angles and for the VK 
(0, 5) and (0, 7) at SZA = 45° for the standard case. Lines with symbols (open squares, SZA = 8°-36°; open 
circles, SZA = 36°- 64°) represent the averaged observed value of the VK (0, 6) band for solar longitude (Ls) 
between 100° and 171° taken from Leblanc et al. (2007). The calculated intensities, when the N2 density is 
reduced by a factor of 3, are also shown. 




Figure 9: The variation of the intensity of the N2 VK (0, 6) emission with respect to solar zenith angle. The 
observed intensity of the VK (0, 6) band is taken from Figure 2 of Leblanc et al. (2007). The calculated intensity 
is averaged- value between 120 and 170 km for the standard case with N2 density in the atmosphere reduced by 
a factor of 3. 
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Figure 10: Intensity variation of the N2 VK (0, 6) band with respect to solar index F10.7 (W/m 2 /Hz) at Mars 
(scaled from the measured value at Earth). The observed intensity of the VK (0, 6) band is taken from Figure 
3 of Leblanc et al. (2007). The calculated intensity is averaged- value between 120 and 170 km for the standard 
case with N2 density in the atmosphere reduced by a factor of 3. 




Figure 11: The calculated limb intensity of the N2 VK (0, 6) band at SZA = 20°. The observed values are 
taken from Leblanc et al. (2007). The calculated intensities are shown for the model atmospheres of Fox (2004) 
(when the density of N2 is reduced by a factor of 2.5) and Krasnopolsky (2002) (the N2 density reduced by a 
factor of 2.1). The intensity calculated by using the electron impact cross sections of Johnson et al. (2005) is 
shown when the N2 density is reduced by a factor of 2. 
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